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Abstract 

We provide an F{R) gravity description of a ACDM bouncing model, without the 
need for matter fluids or for cosmological constant. As we explicitly demonstrate, the 
two cosmological eras that constitute the ACDM bouncing model, can be generated 
by F{R) gravity which can lead to accelerating cosmologies. The resulting F(R) 
gravity has Einstein frame inflationary properties that have concordance to the latest 
Planck observational data. Both the F{R) gravity stability properties are thoroughly 
investigated and also, the gravitational particle production, a feature necessary for 
the viability of the ACDM bounce scenario, is also addressed. As we will show, 
the ACDM bounce model can be successfully described by pure F{R) gravity, with 
appealing phenomenological attributes, which we extensively discuss. 
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Introduction 

Recent observational data regarding the Cosmic Microwave Background radiation mm 
indicate that the scalar perturbations in the early universe were nearly scale invariant, 
and therefore scale invariance is rendered a fundamental requirement that every cosmo¬ 
logical model has to predict to some extent, in order to be considered viable. In most 
cases, scale invariance is achieved by using scalar fields, the perturbations of which [3] are 
scale invariant. Still, no scalar fields have been observed in nature, apart from the Higgs 
probably, therefore a scenario that avoids scalar fields to some extent can be considered 
appealing from a physical point of view. One such cosmological scenario was developed 
by Cai and Ewing [1], in which case scale invariance is achieved by using only ordinary 
cold dark matter, radiation and a positive cosmological constant. Particularly, the model 
describes a bouncing Universe, in which case the Big Bang is replaced by a bounce and 
therefore the Universe is free from the initial singularity. Bounce cosmologies are known 
to be alternative scenarios to the standard inflationary cosmology [3]. For an important 
stream of reviews and papers on bouncing cosmologies, see [5115] . The model developed 
in [1], to which we shall refer to as the ACDM bounce model, makes the important as¬ 
sumption that the equation of state, which describes the perfect matter fluids, changes 
discontinuously and as a consequence of this, the dynamical evolution of the Universe is 
divided in two cosmological eras. The first era is dynamically governed by radiation with 
the assumption that quantum gravity effects are taken into account during this era. Par¬ 
ticularly, the dynamical evolution is described by loop quantum cosmology (LQC) [9HTT]. 
a feature that significantly changes the standard cosmological evolution. The second era 
is governed by cold dark matter and a positive cosmological constant with the change 
between the two eras being continuous, with regards to the scale factor. However, it is 
expected that quantum gravity effects should play an important if not defining role at that 
universe epoch. Hence, the natural question which appears in relation with above scenario 
is: can it be realized within some effective gravity model? It is naturally to expect that 
if such universe can be realized within effective gravity then matter should play minor (if 
any ) role in its occurrence. 

It is one of our main purposes in this paper to provide a pure F(R) gravity description 
of the ACDM bounce scenario, with pure indicating that no matter fields are going to be 
used. In addition, we aim to study the stability of the F{R) gravity we shall reconstruct. 
The F{R) gravity are known to provide consistent theoretical descriptions for cosmolog¬ 
ical scenarios, which ordinary Einstein-Hilbert gravity fails to describe. For informative 
reviews on this vast research stream see m and for important papers consult Refs. [12fB^ 
and references therein. 

It is worth to mention that one of the successes of F{R) theories is the consistent 
description of the late-time acceleration era, with the dark energy finding an appealing 
and self-consistent geometrical explanation. For alternative theories to F{R) gravity that 
also provide a theoretically consistent description of dark energy, see for example [23II28] . 
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Moreover, it is known that modified gravity may successfully realize inflation consistent 
with Planck data (see Refs. [331134] i or even the unihcation of inflation with Dark Energy 
epoch (see Ref. [20] and the reviews |I2|). 

In order to reconcile which F{R) gravity govern the ACDM bounce, we shall make use 
of two quite well known reconstruction techniques developed in Refs. m and m- For 
alternative reconstruction techniques to the ones we shall use, see for example [181119] . In 
addition, for similar studies for other bounces see [22112911^ . The reason we shall use to 
different reconstruction techniques is mainly traced in the particular form of the radiation 
era of the ACDM bounce. For this case we shall use the technique of m which yields more 
reliable results in comparison to the technique of m, as we explain in detail in a future 
work. With regards to the ACDM era described by matter and cosmological constant, we 
shall use the reconstruction technique developed in [15], since it yields exact analytical 
results, without any approximations. In the case of the radiation ACDM bounce era, 
it is not possible to find an explicit analytic solution, so we investigate this case in the 
large curvature limit, most relevant to the LQC era which is governed by large curvature 
quantum operators. As we shall demonstrate, both the eras can be described by F{R) 
gravity that generate accelerating cosmologies. We also study the stability of our solutions 
and we investigate in which case instabilities can occur for our solutions. Interestingly 
enough, one of the two eras is described by an F(R) gravity which when studied in the 
Einstein frame, can be compatible with the latest Planck data. We study in detail the 
cases in which this compatibility can be achieved. In addition to these, we give a brief 
account of the holonomy corrected F(R) gravity [31] corresponding to the matter era of 
the ACDM bounce and we study how the dynamical evolution of the Universe is described 
in this context. Finally, we address the issue of gravitational particle production, which 
is a feature that can render the bounce asymmetric, with the latter being favored by 
observational data. 

This paper is organized as follows: In section 1, after providing a brief description of 
the ACDM bounce and the two eras that it consists of, we make use of the reconstruction 
techniques to investigate which F{R) gravity generate such a cosmological evolution. In 
addition, we study the Einstein frame properties of the F{R) gravity that corresponds 
to the matter era of the ACDM bounce. For the same F{R) gravity we briefly study it’s 
holonomy corrected form in the Einstein frame. The stability of our solutions is thoroughly 
examined in section 2, while the gravitational particle production issue is addressed in 
section 3. The conclusions along with a brief discussion on our resulting F{R) gravity 
picture follow at the end of the paper. 

1 ACDM bounce scenario from F{R) Gravity 

1.1 A Brief F{R) Gravity Review-Conventions 

In order to make the article self contained we review in brief the essential features of F{R) 
gravity in the Jordan frame, and also describe the geometrical background that we are 
working on. There are two approaches in F{R) gravity, namely the Palatini formalism [T2] 
and the metric formalism and we shall work in the context of the latter. The spacetime 
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manifold is assumed to be a pseudo-Riemannian one, which is locally described by a 
Lorentz metric, and specifically the Friedmann-Robertson-Walker metric in the case at 
hand. In addition, the metric compatible affine connection is the Levi-Civita connection, 
which is torsion-less and symmetric. With this connection, the corresponding Christoffel 
symbols are, 

= ^9^^{dti9Xu + d^gxi, - dxg^,u), ( 1 ) 

and additionally, the Ricci scalar is equal to, 

R = ( 2 ) 

The F{R) theories are described by the following four dimensional action in the Jordan 
frame, 

5 = ^ / d'^Xx/^F{R) + Smigfiu, ^m), (3) 

with = SttG and also Sm encompassing all the matter fields present. The metric 
formalism of F(R) gravity is actually materialized if the metric tensor it self is considered 
as the main variable, and by varying action ([3]) with respect to the metric tensor g^^, we 
acquire the following equations of motion 

F'{R)R^,{g) - \F{R)g^, - V^V,F'{R) + gf.,nF'{R) = (4) 

The prime in Eq. ([1]) denotes differentiation with respect to the argument of the differ¬ 
entiated function, that is, F'{R) = dF{R)/dR and additionally denotes the energy- 
momentum tensor of the matter fields. Finally, as we already mentioned, the metric will 
be assumed to be a flat Friedmann-Lemaitre-Robertson-Walker (FRW hereafter), with the 
corresponding line element being equal to, 

ds^ = —df^-|-a^(t) ^ dxf, (5) 

i 

and a{t) being the scale factor. The Ricci scalar corresponding to this line element is, 

R = Q{2H^ + H), (6) 

with H{t) denoting the Hubble parameter H{t) = a/a, and the dot indicates time differ¬ 
entiation. 

1.2 An overview of the ACDM bounce scenario 

The focus in this article is on the cosmological scenario described in detail in Ref. a by 
Cai and Ewing. The scenario itself is quite physically appealing and we shall provide a 
pure F{R) gravity description of it. Before we proceed to the F{R) gravity reconstruction 
of the cosmological scenario, it is worth to provide a detailed description of it, since it is 
of importance to understand the new insights that the F{R) gravity brings along. For 
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details on this scenario, the reader is referred to Ref. 1^. In addition, for similar models 
with distinct cosmological eras, see [32]. 

The cosmological model described in |1|, described the dynamics of a flat FRW cosmol¬ 
ogy with a positive cosmological constant and also with radiation and cold dark matter 
(CDM hereafter) present. The cosmological evolution was divided in two eras, namely one 
described by an effective LQC Hamiltonian, in which the minisuperspace wave function 
is a sharply peaked state, very adequately described by the effective equations of state of 
LQC, and the other evolution era is described by a cosmological constant plus CDM. 

With regards to the effective LQC cosmological era, it was assumed that radiation 
dominates this regime, which is considered to be the high curvature regime, in which case 
the effective equations of motion are given by HEHmEI], 





(7) 


with H the Hubble rate, p the radiation matter fluid energy-density and pc ~ Ppi, the criti¬ 
cal energy-density which is of the order of the Planck energy-density. In the radiation-LQC 
era, the quantum gravity effects are assumed to control the cosmological evolution, as is 
obvious from the holonomy corrected FRW equation ([Tj). With regards to the cosmologi¬ 
cal constant plus CDM era, which in the rest of the paper we shall refer to as matter-CC 
phase, the cosmological evolution is governed by the cosmological constant and CDM. 

In addition, and more importantly, it is assumed that the Universe’s evolution is broken 
into the aforementioned cosmological eras, and that there is a discontinuous transition 
between these two eras, a fact that will be clearly depicted in the equations of state 
describing the two eras. As we explicitly demonstrate, this equation of state discontinuity 
will also be materialized in our pure F{R) gravity description, but we do not include any 
matter fluids in our description. Particularly, as we shall show, the two F{R) gravity 
that correspond to the two distinct cosmological eras, are mathematically distinct, but 
interestingly enough, both produce accelerating cosmologies in the large curvature limit. 
Specifically, the LQC radiation era is described by an R -|- A gravity, and the ACDM era 
is described by an R^ gravity, in the large curvature regime. 

Having assumed a discontinuous evolution of the Universe, in the aforementioned eras, 
let us briefly describe these two eras, the scale factor of which we shall extensively use in 
the following sections. The LQC radiation era is governed by radiation, with an equation 
of state, 

® 

and in conjunction with Eq. dZj) we obtain the LQC radiation era scale factor g]. 


In the earlier epoch, the curvature is lower and quantum gravity effects no longer govern 
the evolution of the Universe, a process which now is governed by the CDM energy-density 
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and cosmological constant, which we denote pcDM and A respectively. In this case, the 
FRW equations are given by the following expression, 

{pCDM + . ( 10 ) 

Denoting the total energy density ptot, h is obvious that 

Plot = PCDM + PA- (11) 


It is also assumed that the equation of state is of the form Ptot = ojptot, with oj = —5 
and 0 < (5 < 1. The parameter 6 is assumed to vary continuously from the cosmological 
constant epoch to the CDM matter epoch, but it is also assumed that when a specific era 
is considered, this is almost constant [Ij. In this context, the total energy density for a 
specific epoch is given by. 


with poff a constant related to the scale factor at the radiation matter era transition, 
which for simplicity we leave it as Peff- Notice that the value <5 = 1 corresponds to the 
cosmological constant epoch {oj = —1), while the value <5 = 0 corresponds to the CDM 
epoch (w = 0). The scale factor for the fluid with equation of state Ptot = ^Ptot, as a 
function of the cosmological time t and of the parameter <5, is equal to, 


a{t) = A{t — 7 ) 3 (i-' 5 ) j 


(13) 


where we have set A to be equal to. 

The parameter 7 appearing in Eq. (|13|) is again related to the transition time from radia¬ 
tion to matter era, but it’s exact definition will play no important role in our analysis, for 
details see [3]. In the following sections, we shall thoroughly investigate which pure F{R) 
gravity can generate the cosmologies described by Eqs. m and m- Special emphasis 
shall be given in the transition epoch between the two eras, and specifically in the high 
curvature regime, when the LQC era ends and the ACDM era starts (and in particular 
the cosmological constant era). This era is of particular importance with regards to the 
F{R) gravity we shall hnd, since it can describe the inflationary era in the Einstein frame. 
We shall study in detail the resulting picture in a following section. 

1.3 ACDM bounce from F{R) Gravity: The matter-cosmological con¬ 
stant phase 

We start our analysis with the reproduction of the matter-CC phase by a pure F{R) 
gravity. Our strategy is to hnd which pure F{R) gravity can produce a cosmological 
evolution with a scale factor equal to the one of Eq. m- We shall use the reconstruction 
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technique of Ref. m, which makes use of the e-fold number N. This is most appropriate 
for the case at hand, since this technique leads to differential equations which can be solved 
analytically. There is an equally useful reconstruction technique developed in m, which 
makes use of an auxiliary field. This technique yields approximate results to the large and 
small curvatures limits, but since we can have exact analytic results with the technique 
developed in |15] . we use the latter. In a following section we shall have a small discussion 
on the issue of choosing the best technique, but let us mention that the two techniques 
yield the same results, if these are used properly. 

The Hubble rate corresponding to the scale factor (fT3|) is equal to, 

2 

3(l-5)(t-7)’ 

and recall that 0 < 5 < 1. The first FRW equation is written in the following way, 

- 18 (4H{tfH{t) + F''{R) + 3 + H{t)^ F'{R) - = 0, (16) 

with F'{R) = and the Ricci scalar R is given in Eq. ([6]). The e-folding number N 

is related to the scale factor in the following way, 



(15) 


-N 

e = — 
a 


(17) 


and by using this variable N, the first FRW equation can be expressed in terms of the 
e-fold parameter N as follows. 


- 18 -h H^{N){H'f -h H^iN)H''{N)) F''{R) (18) 

+ 3 [H^{N) + H{N)H'{N)) F'{R) - = 0. 

In the relation above, the derivatives are defined with respect to the new variable N, that 
is, H' = dH/dN and H” = d?‘H/dN‘^ and the same convention holds true in the rest of 
this section. Introducing the function G{N) = H‘^{N) in equation p8p . the latter can be 
recast as follows, 


- 9G{N{R)) {4G'{N{R)) + G"{N{R))) F"{R) (19) 

+ (^3G{N) + ^G'(iV(R))) F'{R) - ^ = 0, 

with G'{N) = dG{N)/dN and G''{N) = d^G(A^)/dA^^. A crucial point is the fact that 

R = 3G'{N) + 12G(A^), (20) 


and by using this in conjunction with (11711 . we end up to a differential equation which 
will provide us with the reconstructed F{R) gravity which produces the cosmology (fTSl) . 
Indeed, by writing the Hubble rate as a function of the scale factor, 


H = 


2 3 ( 1-7 

3(1^“ 


( 21 ) 
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and by making use of (ED and also recalling that G{N) = we have, 

G{N) = , ( 22 ) 

where A is equal to ■ Thereby, by using Eqs. (I20l) and (f22]l . the e-fold 

number can be expressed as a function of the Ricci scalar, 


with B 


N = - 


1 


3(1 - 5) 


In 



(23) 


—9^4(1 — 6) + 12A, and so the differential equation of Eq. (|19p takes the form. 


9A^3{l-6){3{l-6)-A) 3A(3(1 - 5) - 2) dF(R) F{R) 


0. (24) 


Notice that we assumed that no matter fluids are present so that the differential equation 
(I24D yields the pure F{R) gravity that generates (fTBIl . In order to simplify the notation of 
the equations to follow, we introduce the parameters ai and 02 , which are defined to be. 


oi = - 


9 ^ 23(1 _j)(3(i _ 4) 

R2 


02 = - 


3^(3(1 -6) -2) 


(25) 


The differential equation (j24p is the homogeneous Euler second order differential equation, 
the solutions of which we denote by fi{R) and f 2 {R)- These are equal to. 


h{R) = RP\ f2{R) = R-P\ 


(26) 


where the parameters pi and p 2 are given below, 

-(oi - oi) + (02 - ai)2 -h 2ai -(oi - oi) - (02 - ai)2 2ai 

« =- 2^1 -• « =- 2^1 -• 

and with oi, 02 given in (I25I1 . Hence, the pure F{R) that generates ACDM bounce solution 
of relation (I13p is the following, 


F{R) = ciRP^ + C2RP^ 


(28) 


with Cl, C 2 free parameters. It is worth providing the exact relation of the variables pi and 
P 2 as functions of 6, which is, 


2 1 
- P2 — -• 

10 , 173 - 786 + 95 ^ ’ 3 5 1 173 - 785 + 95 ^ ’ 

1+3(5 Y (1+3(5)^ 2 1+3(5 2 Y (1+3(5)^ 


(29) 


The resulting F{R) gravity varies as a function of the parameter 6 that determines the 
equation of state during the matter-CC phase of the ACDM bounce. We shall investigate 
the form of the F{R) gravity for the limiting values of 6 and also examine the large 
curvature limit which is very essential with regards to the other evolution era of the 
scenario we study, namely the LQC-radiation era. 























The most interesting case from a physical point of view is when —)• 1, which case 
corresponds to an equation of state that describes a cosmological constant. For <5 —)• 1, 
the parameters pi and p 2 behave as follows, 

Pi = 2, lim p 2 —oo, (30) 

< 5^1 

and therefore for 5 —?• 1 and in the large curvature regime, the F{R) gravity behaves as, 

F{R) ~ ciR^, (31) 

since the term RP^ is almost zero as R tends to infinity. We have therefore ended up 
to a very physical appealing scenario in which the large curvature F[R) gravity which 
describes the matter-CC phase is an R? gravity, when <5 = 1. The R? gravity is known to 
have quite interesting inflation properties |33ll34j . This result is of great importance and 
we will discuss it in a later section, when we also have the F[R) gravity which generates 
the LQC-radiation era. From the form of the parameters pi and p 2 as functions of <5, it 
is obvious that in the large curvature regime only the pi term dominates, so the F{R) 
gravity for a general value of <5 in the large curvature regime is, 

F{R)^ciRP\ (32) 

For (5 = 0, which corresponds to the pure matter domination era, the parameter pi is equal 
to, 

Pi = ^(7 + ^ 73 ), (33) 

for which case, nothing interesting occurs. However, having an R? gravity describing the 
large curvature regime, it is worth investigating if there are values of <5 for which we can 
have some overlap with the Planck observational data [2]. In the next section we shall 
investigate the properties of the large curvature F{R) gravity (|32p in the Einstein frame 
and as we explicitly demonstrate, there are values of 6 for which we can achieve exact 
compatibility with the latest Planck data [2]. 


1.3.1 Einstein frame inflation stndy of the ACDM bounce F{R) gravity 


As we evince in this section, the F{R) gravity of Eq. (|28p in the large curvature limit, 
when conformally transformed in the Einstein frame, can yield results that overlap to a 
great extent with the Planck data. To start with, the F{R) gravity (|28p . in the large 
R limit gets simplified and is given by (132 p . A detailed analysis on the Einstein frame 
inflation properties of F{R) gravity was performed in [34] and we adopt the notation of 
this reference. In order to make contact with [M], we identify our parameter pi with the 
following expression. 


Pi 


n + 2 
n + 1 ’ 


(34) 


and also we set ci to be equal to, 


Cl = 7 


1 


n + 1 


n + 2J \4(n + 2) 


l/(n+l) 


(35) 
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since ci is a free variable. Then, the F{R) gravity in the large curvature limit reads, 


F{R) ~ 7 


n + 1 


n + 2J \4(n + 2) 


l/(n+l) 


ti+2 

n+1 


(36) 


Before continuing and making contact with observational data we need to discuss the 
theoretical framework of the Einstein frame. As we discussed in the previous sections, the 
F{R) gravity (|36l) generates the matter-CC phase FRW cosmology with the scale factor 
of Eq. (I13p . so this does not generate inflation in the Jordan frame. Our aim is to study 
the Einstein frame inflation properties of (|36l) . so we assume that in the Jordan frame, 
the metric is an appropriately chosen one, so that when conformally transformed in the 
Einstein frame, it yields an inflation generating de Sitter or at least quasi-de Sitter metric. 
Then by conformally transforming the Jordan frame theory using standard techniques 
], we obtain the following scalar potential, 


Via) = - 


7 (n + 2) 

P 


1 — e 


The slow-roll parameters are defined as. 


1 


e = 


E'(ct) 


2P V V{a) P V V{^) 

during the slow-roll inflation, and for the potential (l371i . these are equal to, 

2 




1 (V'ia) 


(37) 


(38) 


ln—{n + 2)e V 3 i 


/ I r,\ nJ\kP'a , „\ (n 

[n + 2)e V 3 — [n + 2)e 


(39) 


- 1 


V = 


n? + in + 2)e 


(n 


^ 1 + (n + 2)e^” 


rixFik'^c 

3 — (n J- 2je V 3 


These slow parameters at the limit a —)■ —oo, become approximately equal to, 

in? 

e ~ —, r? ~ 2—. 

3 3 


(40) 


So eventually the primordial power spectrum and the spectral indexes of inflation are 
equal to. 


iNn^ ^2 .„2 


n n 

p —-, n, ~ 1 — 2 —, r ~ 16—, 

^ 87r2n2 ’ ® 3 ’ 3 ’ 


(41) 


where we omitted corrections of the order e 3 . The latest observational data constrain 
the spectral indexes as follows. 


= 0.9603 ± 0.0073, r < 0.11, 


(42) 
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so the spectral indexes of Eq. (HU) are consistent with the observational data when n ~ 0.2 
and n ~ 0.1. Recalling relation (I34p and since the parameter pi is related to 5 as follows, 


2 


10 I / 73-78(5+952 ’ 
1+35 V (1+35)2 


(43) 


the values n ~ 0.2 and re ~ 0.1 correspond to 5 ~ 0.95 and 5 ~ 0.97. Of course these 
values are allowed, since 0 < (5 < 1. In the following we shall take into account these two 
values of 5 and the corresponding in each case F{R) gravity, which for 6 = 0.95 reads, 


F{R) ~ ciR 


2-0.166667 


(44) 


while for 5 = 0.97 is equal to. 


F(R) c. 


(45) 


1.3.2 Holonomy Corrected F{R) gravity and the ACDM bounce: An alterna¬ 
tive approach to inflation 

As we already mentioned, the F{R) gravity that describes the large curvature regime of the 
matter-CC phase is of the form ~ RP^ , and during this regime, the quantum gravity effects 
are disregarded. According to scenario of Cai and Ewing, the quantum gravity effects 
effectively modify the LQC-radiation phase. Before proceeding to the F{R) description 
of this cosmological phase, we shall study in this section how holonomy corrections may 
modify the cosmological evolution of the F{R) gravity given in Eq. (|32p . As we shall 
show, the holonomy corrected F{R) gravity in the Einstein frame can successfully describe 
a bouncing cosmology itself, without the need for another cosmological era. We intend to 
give a brief qualitative analysis of the holonomy corrected ~ RP^ gravity in the Einstein 
frame, but we intend to address the problem more concretely and in the Jordan frame in 
a future publication. 

The holonomy corrected F[R) gravity in the Einstein frame were studied in |31] . 
Introducing the Einstein frame holonomy corrections, the FRW equation reads, 

(46, 


with Pc the Einstein frame critical density. Equation (|46p describes an ellipse in the (77, p) 
plane, and the Universe’s evolution is quite simple to describe. Particularly, the universe is 
moving clockwise starting from a contracting phase and proceeding to an expanding phase. 
Both phases begin and end at the same critical point (0,0) and the system bounces off 
only once at the point (0, p). Using the scalar potential in the Einstein frame (1371) . the 
Einstein frame scalar field evolution is equal to. 


a + 377(t + 


dV{a) 

da 


= 0 , 


(47) 
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Performing the transformation = Incr, the Einstein frame evolution of the new 

scalar held V’ is governed by the following differential equation, 


+ (48) 

+ V6(pi^- l)fc ~ - 2{pi - 

with the parameter ”gf” set equal to, 


9 = 


2pi-3 
Pi - 1 


(49) 


The bounce produced from the dynamical system (14811 is symmetric, a fact that makes 
this scenario less appealing, as we shall discuss later on. It is easy to see how a symmetric 
bounce is generated without solving the differential equation (I48|) explicitly, since it pos¬ 
sesses some symmetries that make the qualitative analysis quite easy. Indeed, the orbits of 
the dynamical system (I48p . depicted in phase space by ('ll), ip), are symmetric around the 
Ip, owing to the fact that Eq. (1480 remains invariant under the following transformations. 


t —^ — i, H —y —H 


(50) 


Therefore, the contracting phase {H < 0) orbit {ip{t),ip{t)) under the transformation (I50p 
is transformed to the expanding phase {H > 0) orbit {ip{—t),ip{—t)) and thus, the bounce 
is symmetric. In addition, the Einstein frame energy density p is equal to. 


ip‘^ 1 

^ 2A;2y2 


{{P1A4) ip"^^^ - ^4(^1 A4) i-^>i ip'^'^ 


(51) 


The Hubble parameter H is directly related to the energy-density p by the holonomy 
corrected ERW equation (I46|) and as it is obvious, the Hubble parameter vanishes at the 
point {ip. Ip) = (5,0) and in addition at the curve p = pc- A simple qualitative analysis of 
this result may be done easily and it goes as follows: The universes evolution begins in the 
contracting phase with H < 0 and oscillates around the point (5,0), where the oscillations 
amplitude increases up to the point it reaches the curve p = pc, Sii which H = 0. At this 
point the universe bounces off and it enters the expanding phase, during which H > 0. 
The expansion continues in an oscillating way until the critical point (d, 0) is reached. 

Before we close this qualitative description we have to note two things. Eirstly, the 
p = Pc curve is not simple compared to the gravity curve studied in [3T] and secondly 
the bounce predicted by this scenario is symmetric. A symmetric bounce however is ruled 
out, as was also pointed out by Cai and Ewing [4], so the only way to achieve an asymmetric 
bounce within this framework is only through particle production enhancement. We shall 
address this issue in a later section. 
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1.4 ACDM bounce from F{R) gravity in the large curvature limit: The 
radiation phase 

The Hubble parameter corresponding to the scale factor Q which describes the LQC- 
radiation phase, is equal to, 

327vGpct 

^ 647rGpci2 I o' 

This Hubble parameter will be our starting point for the reconstruction of the LQC- 
radiation phase cosmology. At this point we have to make a crucial remark, related to 
the issue of choosing the most optimal reconstruction method. Particularly, since we are 
interested in the large curvature regime of the LQC-radiation phase, this corresponds to 
early times which means small values of the cosmological time. Particularly, as also noted 
by the authors of Ref. [1], the fact that for small cosmological times the scale factor is 
non-zero is a manifestation of quantum gravity effects. Taking the limit of the scale factor 
([9]) when t tends to zero, we obtain that indeed the scale factor is non-zero and equal to, 

a{t) = 1. (53) 

This is a crucial observation, since the reconstruction method we used in the previous 
section, requires some functional dependence of the scale factor with respect to time, 
so the most optimal method for reconstructing the cosmology described by the Hubble 
rate (l52]) . is the reconstruction method developed in m, which we use in this section. 
Disregarding any contribution form matter fluids, the first FRW equation reads, 

- 18 (4H{tfH{t) + H{t)H{t)^ F”{R) + 3 + H{t)^ F'{R) - = 0. (54) 

The reconstruction method of m, makes use of an auxiliary field cj), so that the action 
© which describes the pure F{R) gravity, is modihed in the following way, 

5 = y d^xV^ {P{(I>)R + Q{4>) + Fmat) ■ (55) 

Practically, the final form of the reconstructed F (R) gravity will be given by the functions 
P{4>) and Q{4>), so the focus is to find these solutions explicitly. The absence of a kinetic 
term for the scalar field in relation (|55p . renders it an auxiliary time dependent degree of 
freedom, so upon variation with respect to cf), we obtain, 

P'{cP)R + Q'i^) = 0 (56) 

where the prime denotes differentiation with respect to (j). Solving this algebraic relation 
with respect to ()>, will yield the function <p{R) and the F{R) gravity can be found easily 
by substituting (piR) to action ([55]), so that F(R) is equal to, 

F{cPiR)) = P{cP{R))R + QicPiR)). (57) 

Therefore, it is of critical importance to modify the FRW equation in such a way so that we 
obtain a differential equation that will yield P{4>) and Q{(p). This differential equation can 
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be easily obtained by varying equation (1551) with respect to the metric tensor. Assuming 
a flat FRW metric, we obtain, 


- 6H^P{^it)) - Qi^it)) - 6H 


dPim 


dt 


= 0 


(58) 


4H + j P(0(t)) + Qi<Pit)) + 2 

By eliminating the function Q{4>{t)) we obtain, 


d^pm)) , dpm)) 


dP 


+ 


dt 


= 0 . 


dt^ 


- 2H{t)P{<j){t)) + = 0. 


(59) 


As it is explicitly proven in the scalar field (p is considered to be equivalent to the 
cosmological time t, so in the following we assume that p = t (see appendix of HZ!)- 
Assuming that the scale factor takes the form. 


a = 


(60) 


with oo being a constant, the differential equation (f5U|l can be recast in the following way, 

(61) 


+ i</'WPm)) = 0 . 


dt^ 


dt 


This differential equation yields P{(p), and by using the resulting P{(p), we may get the 

Q{^), 


QiP) = -6g'iPfPiP)-6g'{P) 


dPjP) 

dp 


The Hubble rate (I52p can be written in the following form, 

h{t) 


H{t) = 


t 


with h{t) being equal to. 


h{t) = 


hjqt^ 


l + qP' 

and where we defined hj and q to be equal to, 

, 1 32ttGpc 

hr = q = - 

■'2 3 


(62) 


(63) 


(64) 


(65) 


The function h{t) appearing in Eq. (1641) is a slowly varying function of time, a property 
that will significantly simplify the problem at hand. Indeed, the function h{t) V z G M 
satisfies the following relation, 

( 66 ) 


lim 


t-^oo h{t) 
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We therefore assume a particular form for the function appearing in relation (I60p which 
is, 

5f((/)) =/i((/.) In , (67) 

with 00 some integration constant. Since the function /i(0) is slowly varying, we can ignore 
it’s derivatives in the following calculations. Using the functional form of g'(0), the Hubble 
rate is equal to, 

^+ , ( 68 ) 

and owing to the fact that the function h{t) is slowly varying, the derivative h'(t) can be 
safely ignored, and thereby the Hubble rate simplifies to. 





(69) 


which is exactly equal to the LQC-radiation phase Hubble rate given in Eq. (|6.SD . This is 
a crucial point in our analysis, and the validity of the method is verified by exactly this 
coincidence. By using Eq. (l671l and ignoring the derivatives h'{t),h"{t) the differential 
equation (I6ip can be cast in the following form. 


d^P(0(t)) /i(0) dP(0(f)) 

dt^ 0 dt 


2h(0) 

02 


P(0(t)) = 0. 


(70) 


In order to find the P(0(P)) function from the above equation, we must first find the 
exact functional dependence of 0 as a function of the Ricci scalar, which can be easily 
done if we calculate the Ricci scalar R by taking into account relations ([6]) and (1631) and 
(I64p . Combining these, the final result reads. 


P(0) ~ 


6 (-/i(0) + 2/i(0)2) 

02 


(71) 


Notice that we ignored the higher derivatives of h{t) in order to obtain Eq. (I7ip . Using 
the exact form of /i(0) given in Eq. (16411 (recall that 0 = t), solving Eq. (j71ll with respect 
to 0^, yields the result, 


_ 2 2 hf 

~ P ~ ^ 3R 


(72) 


/ I/O 12x2^/3 2 21/3^ 2^/^qhf 

V ^ R 3 3 

4 X 2V3g2/^, 8x2^/^qh‘l 2bx2^/^q^h) Ax2^l^q^h?r i/s24 1 

- - — - - - H- - - - + 4 X 2^''^q^h f H- -jr 

P O ' op -TP / I /O 


32V3g2^ 


ao 


+ aiR + 02^2 + aaPS + ^ ^^R^ + p^R^ + p^R^ + ^^R^^ 


1/3- 
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For simplicity, the coefficients Oi, I3i are given in appendix A. Using Eq. (1721) . we can solve 
the differential equation (1701) and express the solutions as functions of the Ricci scalar. 
Then, by using Eq. (j57p we can have the reconstructed F{R) gravity. The general solution 
to the differential equation ([70l) is [15], 

_ , , ^(0) —l+V +6h(4>) + l h(cl>)-l — y/h((li)‘^+6h(4>) + l 

= cKk ^ + C2<i> ^ , (73) 


and by using this we may easily find Q{4') by substituting (IU7|) in equation ([H2I1 . in which 
case we obtain, 


Q{(p) = -Qh{^)ci h{4>) + 


— It T 6 /i((^) T 1 \ ^ _ 2 


(74) 


- 6h{(p)c2 h{(p) T 


h{(f>) — 1 — y^/i((?!))2 T 6h{4>) T 1 \ ^ 


Since we are interested in the large R regime of our approximate method, in the following 
we shall examine the solutions we found in the large R limit only. In fact, it is an exercise 
of academic interest only to examine the small curvature limit, since in the context of 
the cosmological scenario we are studying, this era is governed by the CDM fluid, so it is 
senseless to examine the low curvature limit. 


1.5 F{R) Gravity in the Large R Limit 

Having solution (1721) at hand, we may easily obtain the large R limit, which is, 




1. A. 

3q^ R 


so finally we get, 


IR-Ai 

3qR 


which holds true when i? > ^i, where we have set .Ai to be equal to, 

4 X 2V3g + ^2!!^ _ + 4 X 2F3q^hj + {2 - hf + 2h})al^^ 


Ai — 


1/3 


(75) 


(76) 


(77) 


In addition, the parameter 03 can be found in appendix A. We can find the function 
P{4>{R)) in the small 4> limit (or large R limit equivalently), by recalling that the function 
h{t) has the following limiting value as (j) approaches zero, 


lirn h{t) = 0. (78) 

t^o 
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This observation significantiy simpiifies the caicuiation, since in that hmit P{(f)) reads, 


P((/>) = Cl + C24> ^ 

In virtue of Eq. (ITHI) . the term is equai to, 


A-i 


3qR 


R — 

so eventuaiiy, the function P{4>) of Eq. (f7^ is equai to, 


P{(t)) = C1 + C2 


2>qR 


(79) 


(80) 


(81) 


R — 

In the same vain, the function Q{(p{R)) is found to be approximately equal to zero, 

Qi<P{R)) ~ 0. (82) 

Combining Eqs. (|81|) and (I82|) . the final form of the reconstructed F{R) gravity of Eq. 
(IFTli reads. 


F{R) ~ cii? + C 21 


' 3qR^ 
R — J\.i 


(83) 


We can further simplify the resulting F{R) gravity, by expanding the square root in the 
large R limit, in which case we get. 


F{R) 




which is rewritten. 


F{R) ~ ^ci + C2\/^^ R + 


C2^/^Al 


(84) 


(85) 


We can choose the coefficient of R in the above expression to be equal to one, that is. 

Cl+C2 1 /^= 1 , (86) 

and defining A in the following way, 

C2^/^Al 


A = 

2 

we end up to the final form of the reconstructed F[R) gravity, 

F{R) ~ +A, 


(87) 


( 88 ) 


which is Einstein-Hilbert gravity plus cosmological constant. This is a quite intriguing 
result, since the cosmological constant can describe inflationary dynamics of an expanding 
universe. Let us here recapitulate what we found in this section. We investigated which 
F{R) gravity can generate the large curvature LQC-radiation cosmological scenario, which 
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corresponds to the scale factor (1521) and we ended up to an Einstein-Hilbert gravity plus 
a cosmological constant. Notice that the value of cosmological constant A can be of the 
Planck energy order, since the parameter q is related to the Planck energy-density (see 
Eq. ()65p l. Combining this result with the one we obtained in the ACDM study, which 
resulted to an gravity we have the physically appealing picture in which, when the 
large curvature limit of the two cosmological scenarios is considered, both reconstructed 
F{R) gravity result to gravities that actually can generate inflationary solutions. And 
most importantly, in our case no matter fluids are present, so the inflationary dynamics is 
a result of pure F{R) gravity. We shall further discuss this result in a later section. 


2 Stability of F{R) gravity describing the ACDM bounce 
and the radiation bounce 

Before proceeding to some phenomenological applications of the modified gravity descrip¬ 
tion for the ACDM bounce we provided, it is of critical importance to check the stability 
of our solutions. In this section we address this issue using the formalism developed in 
Ref. [I?]. We start off with the F{R) gravity given in Eq. ([28]) . which describes the ACDM 
phase of the bounce solution. 


2.1 Study of the F{R) stability for the matter-CC phase 

Eor the modified gravity of Eq. (|28p we used the reconstruction method developed in 
for which, the most appropriate stability check method is the one firstly done in 
which we employ in this section. We perform a perturbation of the form. 


G{N)=g{N) + 6g{N), 


(89) 


and we insert this in Eq. (fT^ noticing that g{N) satisfies equation ([T9]). Consequently, 
the perturbation Sg{N) satisfies the following equation. 


giN) 


<FF{R) 


dR2 


R=3g'{N)+12g{N) 


6"g{N) + [3ff(iV) {Ag'{N)+g"{N)) 


R=3g'{N)+12g{N) 

(90) 


+ 1 MN) - i/(iV)) 21^ 


+ 


12g{N) {4g\N)+g"{N)) 


R=3g'{N)+12g{N)i 

d^F{R) 


6'giN) 


+ {-4g{N)+2g'{N)+g"{N)) 


dR3 
d^F{R) 


R=3g'{N)+12g{N) 


dR2 


+ 


1 dF{R) 


R=3g'iN)+12g(N) 3 dR 


R=3g'{N)+12g{N)i 


6g{N) = 0. 


Using this equation, the stability of the solution under small perturbations of G{N) may 
be directly obtained. Indeed, the stability conditions read. 


Ji 


6i4g'{N)+g"{N))F'"{R) g'jN) 

F"{R) ^ g{N) 


(91) 
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and in addition, 


J2 


m{Ag'{N)+g"{N))F"'{R) 6g'{N) 3g"{N) F'{R) 

F"{R) ^ g{N) + g{N) ^ g{N)F"{R) 


(92) 


From Eqs. dHQD, (mi), (mD and in conjunction with Eq. (USD, we can check the stability of 
the solution (|28p . under the small perturbation ()89p . in which case the stability conditions 
become. 


Ji = 6 + 3(1 - 5) - 23(1 - <5) - 4 + 3(1 - 5)) (-2 + pi)(-l + pi)pi 

(93) 

+ 3^^C2( - ^( - 4 + 3(1 - 5))e-3(i-'5)^y'(-2 + p2){-l + P2)p2) 

X 3^1 Cl ( - - 4 + 3(1 - ,5))"*'(-1 + Pi)pi 

+ 3^2c2 ( - - 4 + 3(1 - ,5)) (-1 + P2)p2 > 0, 


and in addition, the parameter J 2 reads, 

J 2 = -12 - 63(1 - 5) + 27(1 - 5f + (^ci - 37l3(l - Pi 

(94) 

+ 02(12246-3(1-5)^ - 3^3(1 - (j)e-3(i-^)5v)"'^"V2)^-'Qr' 

+ 36( - 4243(1 - (5)e-3(i-'5)i^ + 24c^e-3(i-'5)i^) 

X (ci(l2Ae-3(i-'5)i^ - 3^3(1 - (5)e-3(i-'5)i^) {^-2 + pi){-l + pi)pi 

+ C2(l27le-3(i-'5)iv -3^3(1 -(j)e-3(i-^)5^)"^^^'(-2 + p2)(-l + /02)/02) x > 0, 

where Qi is equal to, 

Qi = Cl (l2yle“’('“'’)'' -3/13(1 +pi)()i (95) 

+ cj ( i 2 / 1 c - 3 C -')" - 3 / 13{1 - 5 ) c - 3 C -')")1 + 


and also A = Having the stability parameters Ji and J 2 we can directly check 

the stability conditions for the most interesting values of 5 we found in the previous 
sections. Particularly, for <5 = 1, in which case the F{R) gravity becomes R^, and for 
6 = 0.95, 6 = 0.97, for which values, the corresponding Einstein frame theory produces 
inflationary parameters which have exact concordance with the observational data. For 
(5 = 1, the stability parameters become. 


Ji — 2 > 0, J 2 — 0, 


(96) 
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thus the E? does not satisfy the second stability solution, a fact that shows that this 
solution is unstable under perturbations. However, when 5 takes the value (5 ~ 2 but not 
exactly two, stability is ensured. In the same way, the stability parameters for 5 = 0.97 
become, 

/ \ 

0.803333 


<7i — 3 


2.03667 + 


V 


2.8063x1048c1 ^ 


1 + 


“0 _ 


> 0 , 


c2 


7 


with regards to Ji which is obviously positive. As for J 2 , this reads. 


(97) 


J2 = 


e-0-0-(4.7110--C2(^)^^-7.947c, 

, 41.61n 


a, 


f2.2810-49c2 (2.2810-49 c2 + 0.64C1 

y-O.OQiV ^ 29.2 7e_::^_^V2.67 


„^0.18g0.187V ( _2.1410-49 c2 


155“ 


H-^- 

a[j-99 (^2.2810-49C2 (2.2810-49 c2 


+ 8.4ci 


ao 


15155“ 


,-0.09]V ^ 28.2 /p-0.09]V ^ 41.61\ ’ 

+ 0.64ci ^-:n 5 ^ 


15155“ 


'‘0 / \ “0 
with the most dominant term for any value of the e-folds number N being. 


(98) 


J 2 — 


«o 


0.18„0.18Af 


8.4ci 


e-0 09^ \ 42.61 


ao 


5155" 


aO-99 (2.2810-49C 2 (2.2810-49c 2 + 0.64ci 


41.61 


> 0 


(99) 


which is obviously positive. Therefore the solution for 5 = 0.97 is perfectly stable. The 
same applies for 5 = 0.95, in which case we have, 


/ 


Ji = 3 


\ 


2.03667 + 


0.803333 


V 


2.8063x1048Cl { o.o^ ) 


> 0 , 


C2 


which is positive, as for J 2 the leading order contribution reads. 


( 100 ) 


J 2 — 


8.46aFSeO'48^ci 


\ 42.61 

“^^5U5- 


ag-09 (2.2810-49 c 2 (2.2810-49 cs 


,-0.09JV ^ 28.2 / p-0.09iV ^ 41.61 

+ 0.64ci 


15155“ 


> 0 , 


( 101 ) 


which is also positive. Therefore, as can be easily checked, for all values of 5, except for 5 = 
2, stability is ensured. In addition, from an observational point of view phenomenologically 
most interesting cases S = 0.95 and 5 = 0.97, provide stable F{R) solutions. 
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2.2 Study of the F{R) stability for the radiation bounce case 

We now examine the stability of the LQC-radiation generating F{R) gravity given in Eq. 
(IHSll . Since this F{R) gravity is the result of a reconstruction method, different to the 
previously used, we shall use the stability method developed in m- We start from the 
following equation. 




dt^ dt 

which can be rewritten in the following way, 


( 102 ) 


^^(*1 +s(‘P)^)pw-o, 


(103) 


and after some calculations, it can be cast as follows. 


((^)^ - 0 - KT" ^ - »■ 


(104) 


We define the function 6 to be equal to, 


<5 = 


dt 


(105) 


This parameter 6 represents actually the way that perturbations behave for the solutions 
we found, since it practically measures the deviation of the parameter (j) from t, which we 
considered to be identical. Using the parameter 5, we can recast Ea. (ll04l) in the following 
way, 

^ = -u}{t)6, (106) 

where uj{t) stands for. 


u{t) = 2 


dp 


-g'{^)^+g''{<P)P{^) 

™ + 25'(</.)P((/.) 



(107) 


If w > 0 for the solution P{(p) we found then stability is ensured, since the perturbation 
becomes small. In the contrary case, instability occurs, since the perturbations are large. 
We shall examine the function P{(p) that corresponds to the large curvature limit, since 
this is the most interesting case. Eor the function g{4>) given in relation (1671) . u becomes, 


UJ 


i^) = 


-P{4>) + H4>)h'{(l>)) + 2P{4>) (-^ + ^ + H<t>)h''{4>)) + P''{ct>) 


2P(</,)(^ + ln(0)h'(0)) +P'(<))) 


(108) 


21 




















and by neglecting the higher derivatives of u becomes, 

+2P(</,) +P"(^) 






(109) 


In this case, by using P{(i>) as given in Eq. (|79h . the form of /i((/>) given in (I64p and 
substituting in Eq. (11071) . we obtain. 




2C2 


2qcihf 


9 H 14^02- 




( 110 ) 


q4>cihf 


qc2hf 






2qc2hf 


0(l+9« f-f 


which for small (p (which corresponds to large curvatures) becomes, 

w(</>) -§+ #, (111) 

9 

where q is defined in Eq. (|65ll . Therefore uj is positive and thus stable for qcj)^ > 2. On 
the contrary w is negative for qcjp' < 2 and in this case the solution is unstable. 


3 Gravitational particle production 

As was also pointed out in Ref. [1], asymmetry of the bounce is required in order these 
cosmological bounce to be viable. One feature that can cause asymmetry of the bounce and 
that guarantees viability of the model, is particle production during the bounce. In this 
section we shall examine the gravitational particle production issue, adopting the research 
line of Refs. [35]. As is well known |35|, curvature oscillations can generate gravitational 
particle production. What mainly interests us is to examine the gravitational particle 
production for the matter-CC phase, which as we found is described by the F{R) gravity 
given in Eq. (1281) . since the LQC-radiation phase results to an P + A gravity with known 
gravitational particle production properties [36] . 

In the large curvature limit, the F{R) gravity that generates the ACDM phase is 
an p2 gravity, in which case there is sufficient particle production [35] to guarantee the 
asymmetry of the bounce. We shall not go into details for this F[R) gravity, since this 
issue was addressed in full detail in [35], but we shall be interested in the small curvature 
limit of the ACDM phase, in which case the F{R) gravity is approximately equal to, 

F{R) ~ C2RP^ (112) 
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since p 2 is a negative number. Recalling the details with regards to the values that p 2 
takes, since 0 < <5 < 1, the parameter p 2 takes the following values, 


-oo < p 2 <-0.128 (113) 

where the limiting values are obtained in the following limits. 


lim p 2 = —0.128, lim p 2 = —oo (114) 

5^0 < 5^1 

Let us briefly recall here how the ACDM bounce scenario works. At the beginning we 
have the bounce and when the curvature is quite large the radiation phase governs the 
dynamical evolution of the Universe. After this phase, the cosmological constant and cold 
dark matter governs the dynamics. At first, and when the curvature is too large, the 
cosmological constant phase occurs, which corresponds to (5 = 1, and as the curvature 
lowers, the matter fluid governs the expansion. In the latter case, 5 approaches zero. 
We shall study the gravitational production during that phase and for particular limits. 
Notice that in the end of the matter fluid matter-CC phase, a new bounce occurs, so this 
is a late time era, with time considered to start when the previous bounce had occurred 
in the past. Following [35], the equation that will reveal the particle production rate is 
the following, 

2>UF^r -R + RFr - 2F{R) = 0 (115) 

where we assumed that matter is absent. In Ref. [35|, it was assumed that the spacetime 
is a flat Minkowski one, but here we assume a flat FRW metric of the form ([5|), in which 
case the Ricci scalar depends only on time. In this case, Eq. (I115p is simplified and by 
introducing the new variable y = RP^~^ and also using the form of the F{R) gravity given 
in Eq. (11121) , Eq. (|115p becomes, 

3p2d^y + SHdty - _ 2)yP2/(p2-i) ^ q 

Keeping only the dominant terms in the small curvature limit, Eq. (11171) gets simplified, 

:ip2d^y + SHdty + {p2 - = 0 (117) 


where H is the Hubble rate m- We shall study the above equation in various limits in 
order to see how particle production behaves in the small curvature limit. As was proven 
in [35], the rate of gravitational production is given by. 


Ppp ^ 


11527r 


(118) 


where is the amplitude of the curvature oscillations and w their physical frequency. 
Note that IS.r is a slowly varying function of t in general. The method used in [35| required 
that the solution y is definitely written in the following form, 


R(t) = yV(P2-i)(t) = Afi(t)sin(a;t) 


(119) 
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since a general analytic solution can be quite difficult. In our case we shall find analytic 
solutions of (|117l) . with one of them being exactly of the form (I119p . without any approx¬ 
imation. The other solution we shall present is analytic but we fit the resulting function 
in such a way that we end up to a function of the form ()119l) . We start our investigation 
with the assumption that (5 —)• 1 and also for times such that, 


2p2 

1-6 


> t 


( 120 ) 


Practically, this means that we are dealing with times after the LQC-radiation phase, and 
during the start of the ACDM phase. In this case, by keeping the dominant terms in Eq. 
(ESI), the latter becomes, 

3Hdty + {p 2 - 2 )/ 2 /(p 2 -i) = 0 (121) 


which has a simple solution, 

-(t-ci)(l-5) 

y{t) = ce -^- (122) 

The function (jl22p if written in the form (I119p . which can be done by a numerical fit, 
yields the following rate of gravitational particle production. 


Ppp ~ 6079 


(123) 


where we used <5 = 0.99 for which case p 2 — —33.19. Notice that this rate strongly 
depends on the fraction (1 — 6 )/p 2 and on the arbitrary parameter c. We proceed to 
another interesting limiting case, which has an interesting analytic solution. Particularly, 
we are interested in the late time evolution, and particularly for the cosmological time 
being, 

t » ^ (124) 

In this case, and as d —)• 1, the differential equation (|117p takes the following form, 

3p2dly + {p 2 - 2)y^2/(P2-i) = g (125) 


which is the Emden-Fowler differential equation. Setting, 


m = 


P2 


/02 - 1 ' 

the solution of (I125p for m ^ —1 becomes. 


A = - 




1 1 


1 


1 


P2 — 2 

3/92 


2A 




2’1-t-m’ 1 + m’ C'a(l-|-m,)' 


) + Ch 


(126) 


(127) 


with El being the Gauss hypergeometric function and Ca, Ch arbitrary integration con¬ 
stants. When (5 —)• 1, the parameter m is approximately equal to m ~ 1, since only p 2 
dominates. In addition, A = —1/3 therefore, in this limit, the solution (11271) is written, 

t + Ch (128) 
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and by making the replacement z 
the following functional form, 


^y, the Gauss hypergeometric function Fi has 


+ z ^arcsin(z) 

Using (11291) . solution (11281) in terms of y{t) becomes, 


(129) 


t = arcsin(-—y) + Gj, = 


(130) 


so hnally y(p) is exactly equal to. 




(131) 


Note that we arrived at the solution (I13ip using only analytic methods, without any 
approximations. Notice that yip) is exactly of the form (|119l) . so in the present case, the 
amplitude and the frequency of the curvature oscillations are. 




1 _ 1 


(132) 


Consequently, the rate of the gravitational particle production is, 

(3Cj)ll52,r 


(133) 


Since the parameter Ca is a free parameter of the theory, it can be chosen to be quite 
small, so that the particle production rate is as big it is required in order to obtain an 
asymmetric bounce. Therefore, in our F{K) gravity description of the ACDM bounce, 
sufficient gravitational particle production is ensured to guarantee an asymmetric bounce, 
during of course the matter-CC phase we described above. 


Discussion 

We provided a pure F{R) gravity description of the ACDM bounce scenario that was 
developed in Ref. [3j, without the need for perfect matter fluids to govern the Universe’s 
dynamics. The ACDM bounce scenario consists of two distinct cosmological eras, namely 
the radiation dominated and the matter-CC phase, and is therefore based on the disconti¬ 
nuity of the equation of state between the two eras. Using very well known reconstruction 
techniques, we were able to find which pure F{R) gravity can generate each cosmolog¬ 
ical era. In the case of the radiation phase, the curvature is considered to be large, so 
the scale factor is a{t) ~ (aP + 1)^^^, a result that is obtained by using LQC consid¬ 
erations. In the large curvature regime, such a cosmological expansion is generated by 
an F{R) ^ R + A gravity, plus non dominant curvature terms in this approximation. 
In addition, the matter-CC era is generated by a power law F(R) gravity of the form 
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FiR) + C 2 -R^, with Pi being numbers related to the details of the scale factor. 

In the large curvature regime, this F(R) gravity is described by an R^ gravity, while as 
curvature lowers, the F{R) gravity takes the form with pi < 2. In this high curvature 
regime we found two particular values of the parameter pi, for which the F{R) gravity, 
when studied in the Einstein frame, yields results that have concordance with the latest 
Planck data on inflation. 

Interestingly enough, in the context of F{R) gravity we found a solution to a problem 
that the ACDM bounce model of [1] was confronted with. Particularly, in order that 
the ACDM bounce model is considered viable, the bounce must be asymmetric. As was 
pointed out in [3], one way to achieve this is through particle production during the bounce 
and in the framework of F(R) gravity this process occurs naturally. As we explicitly 
demonstrated, particle production is particularly enhanced during the matter-CC phase, 
thus the asymmetry of the bounce can be ensured with the F{R) gravity description of 
the ACDM bounce. Note finally that using the same method one can reconstruct ACDM 
bounce universe in F{G) and F{T) gravity (see Ref. [22]) this will be done elsewhere. 

A very important remark is in order. We have to mention that is of fundamental 
importance to explain why nature should select the discontinuous change in the equation 
of state. This is closely connected to the LQC effects, which effectively modify the large 
curvature era and govern the early time evolution of the bounce. This discontinuity of the 
equation of state appears in the F{R) gravity description we provided, since in the large 
curvature regime, the LQC radiation phase is generated by an almost Einstein-Hilbert 
gravity, while the large curvature regime of the matter-CC phase is governed by an R^. It 
would be quite interesting to find a natural explanation of this discontinuity in the context 
of Jordan frame F{R) gravity, an issue we hope to address in a future publication. 

It is important to discuss the possibility that the F{R) realizations of the ACDM 
bounce we just presented, can be distinguished from other theoretical implements. In the 
case of Ref. |1|, it was explicitly demonstrated that the running of the spectral index is 
negative, and therefore it can be tested from CMB observations. In the case of F{R) 
gravity now, in principle a difference may come from the calculation of the spectral index, 
since it is calculated in a different way. But before getting into the details of this, let 
us mention that prior of distinguishing F{R) gravity prediction, it is important to find 
seeds of a bouncing cosmology to the observational data, since there exist many bouncing 
scenarios. In addition, a quite interesting scenario was presented in Refs. where a 
contracting bouncing phase preceded the slow-roll inflationary phase. With respect to the 
latter scenario, this could have observable effects in the CMB anomalies. By determining 
the exact evolutionary scenario, in the context of F{R) gravity, the observational indices 
can be calculated explicitly in two frames, namely in the Jordan frame |38| and in the 
Einstein frame |3l|. The Einstein frame calculation involves a scalar field on which the 
slow-roll conditions are imposed. Interestingly enough, when the calculation is performed 
in the Jordan frame, it is possible to have results that can directly be fitted to the ob¬ 
servational data |38|, without the need for imposing the slow-roll conditions. Indeed, in 
the Jordan frame, if no matter is present, as in our case, there is no scalar field and 
hence contact with observations can be done by using the technique of maximum likeli¬ 
hood. In this way, the parameters of the theory can be appropriately adjusted, so that the 
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corresponding observational indices are produced. The observational indices have quite 
complicated form to be explicitly presented here, but can be found in the Appendix C of 
Ref. [38] • Notice that pure F{R) gravity models in the Jordan frame are in principle less 
restricted from the corresponding Einstein frame counterparts, owing to the fact that the 
parameters can be appropriately chosen. Hence, in order to distinguish these from other 
theoretical descriptions, other observational quantities must be examined, in addition to 
the aforementioned observables, where the differences could be signihcant, like for example 
in the growth index or similar quantities. Another possibility to distinguish predictions of 
modified gravity from General Relativity, may be related with the description of compact 
massive (neutron) stars but this goes beyond this work. 

Finally, of the most sound results obtained in Re. [I] is the existence of red tilted 
long wavelength perturbation modes. It would be therefore very important to study the 
long wavelength perturbations in the context of Jordan frame F(R) theories. A study 
of perturbations valid to a certain limit was done in [39|. Of course the Einstein frame 
analysis can be dealt with standard techniques, but if someone addresses the full problem 
in the Jordan frame, this study can be quite difficult. We hope to address these issues in 
a future work. 
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APPENDIX A: Detailed presentation of polynomial coeffi¬ 
cients 

In this appendix shall give in detail the polynomial coefficients with i = 0, ..3 and 

j = 0, ...5, which appear in Eq. (j72p . Particularly, the aj’s are the coefficients of the 
following polynomial: 

y9(R) = CHq T CX\R CX 2 R? ~\~ Oi^R? (134) 

while the /3j’s are the coefficients of, 

P{R) = + hR^ + + kR^ (135) 
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These are given in detail below, 


ao = 432g® - 216q^hf + A68q^hj - U6q^hj + 156q^hj - 24q^hj + IGq^hj (136) 
«! = 216g® — ISq^hf — 75q^hj + SOq^h^j — ASq^hj 
02 = 36 ( 7 '^ + Sq'^hf + 30q‘^h‘j 
03 = 2 q ^ 

/32 = 456Mq^°hj - 15228q^°hj + 3l725q^^hj - 5076q^°hj + 5076g^°/i^ 

/Js = 23004:q^h} + bMq^h) - I27>28q^h) - 216q^h} 

13a = 3861q^h} + 540q^h} - W8q^hj 
135 = 2lQq'^hj 
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